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ABSTRACT 

Context. X-ray spectra of many accreting pulsars exhibit significant variations as a function of flux and thus of mass accretion rate. 
In some of these pulsars, the centroid energy of the cyclotron line(s), which characterizes the magnetic field strength at the site of the 
X-ray emission, has been found to vary systematically with flux. 

Aims. GX 304-1 is a recently established cyclotron line source with a line energy around 50 keV. Since 2009, the pulsar shows 
regular outbursts with the peak flux exceeding one Crab. We analyze the INTEGRAL observations of the source during its outburst in 
January-February 2012. 

Methods. The observations covered almost the entire outburst, allowing us to measure the source’s broad-band X-ray spectrum at 
different flux levels. We report on the variations in the spectral parameters with luminosity and focus on the variations in the cyclotron 
line. 

Results. The centroid energy of the line is found to be positively correlated with the luminosity. We interpret this result as a manifes- 
tation of the local sub-Eddington (sub-critical) accretion regime operating in the source. 
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1. Introduction 


In accreting binary pulsars, matter from the normal stellar com- 
panion is transferred to a highly magnetized (B >10 12 G) neutron 
star. In the vicinity of the accretor, the gas flow is disrupted by 
the neutron star’s magnetic field and channeled towards the mag- 
netic poles, where most of the X-rays originate. The physics and 
the structure of the X-ray emitting region(s) above the neutron 
star surface are still highly debated (see, e.g., Becker & Wolff 
2007; Farinellietal. 2012). Since the matter hitting the accre- 
tor’s surface is highly ionized, the magnetic field strength is a 
crucial parameter determining the physical processes inside the 
emitting region and the formation of the observed X-ray spec- 
trum. A direct way to assess the B-field strength at the site of X- 
ray emission is the measurement of the cyclotron resonant scat- 
tering features (CRSF or cyclotron lines) in the X-ray spectrum 
of a pulsar. These features appear as absorption lines, caused by 
the resonant scattering of photons off the electrons in Landau 
levels (see, e.g., Triimper et al. 1978; Isenberg et al. 1998). The 
energy of the fundamental line and the spacing between the har- 
monics are directly proportional to the field strength. 


In some accreting pulsars, the energy of the cyclotron line 
has been found to vary with luminosity, apparently due to 
a displacement of the line formation region. Such variations 
of the line energy have been reported for V 0332+53 (e.g., 
Tsygankov etal. 2010), 4U 01 15+63 (e.g., Tsygankovetal. 
2010), HerX-1 (Staubert et al. 2007; Klochkov et al. 2011), and 
A 0535+26 (Klochkov et al. 201 1). The luminosity-dependence 
of the cyclotron feature has strong implications for the physics 
of the X-ray emitting region as discussed, e.g., by Staubert et al. 
(2007), Klochkov et al. (2011), and Becker & Wolff (2007). 

GX 304-1 is a recently established cyclotron line source 
(Yamamoto etal. 2011). It was discovered in a balloon exper- 
iment in 1967 and subsequently identified as an X-ray pulsar 
with a period of ~272s (McClintocketal. 1977). The system 
contains a Be-type optical companion (Mason et a!. 1978) and is 
located at a distance of ~2.4kpc(Parkes etal. 1980). Since 1980, 
GX304-1 has remained in a quiescent state, showing no out- 
bursts. Starting from 2008, when the source was detected with 
INTEGRAL (Manousakis et al . 2008) , GX 304- 1 “resumed” its 
activity exhibiting outbursts with a period of -132 5 d. 

The energy of the cyclotron line in GX 304- 1 was measured 
with Suzaku and RXTE to be around 52keV by Yamamoto et al. 
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Thble 1. INTEGRAL observations of the GX 304-1 outburst in 
January-February 2012. 


Rev. 

Obs.ID 

Mid. MJD 

Exposure [ksec] 
JEM-X IBIS SPI 

1131 

09400230006 

55944.0 

6T6 

427 

68.6 

1132 

09400230007 

55947.0 

42.4 

31.9 

36.6 

1133 

09400230008 

55950.0 

- 


10.7 

1134 

09400230009 

55952.8 

73 

25.4 

37.8 

1135 

09400230010 

55955.7 

- 

6.7 

25.1 

1136 

0940023001 1 

55958.7 

36.9 

28.1 

32.9 

1137 

09400230012 

55962.0 

78.1 

59.7 

78.4 

1138 

09400230013 

55965.0 

60.7 

45.2 

523 


(201 1). These authors demonstrate that the data taken at different 
flux levels show an indication of a positive correlation between 
the line energy and the X-ray flux, although at a low signifi- 
cance level. This made GX 304—1 a good target for a luminosity- 
dependent study of the cyclotron feature. In this work, we 
present the analysis of INTEGRAL data (see next section for 
the mission description) taken during the January 2012 outburst 
of the source. The data reveal a positive correlation between the 
cyclotron line energy and the source flux, as well as the variation 
in other spectral parameters during the outburst. We discuss our 
finding in the context of a model assuming that different accre- 
tion regimes can operate in a particular pulsar depending on its 
X-ray luminosity. 

2. Observations and data reduction 

At the beginning of January 2012, GX 304- 1 entered an outburst 
as reported by Yamamoto et al. (2012). The outburst was moni- 
tored by the International Gamma-Ray Astrophysics Laboratory 
- INTEGRAL (Winkler et al. 2003), starting at MJD ~55943.5, 
when the source flux in the -20-80 keV range was ~250mCrab, 
through the maximum of the outburst, when the source flux ex- 
ceeded one Crab, to MJD -55965.5, when the flux dropped 
to -100 mCrab. INTEGRAL performed one observation every 
satellite orbit (about three days), with a typical exposure of a 
few tens of kiloseconds each. In total, eight observations were 
performed. The INTEGRAL scientific payload contains three 
X-ray instruments: (i) the imager IBIS sensitive from -20 keV 
to a few MeV (Ubertini et al . 2003; Lebrun et al . 2003); (ii) the 
spectrometer SPI sensitive in roughly the same energy range as 
IBIS (Vedrenne et al. 2003); and (iii) the X-ray monitor JEM-X 
operating between -3 and -35 keV (Lund etal. 2003). Table 1 
summarizes the INTEGRAL observations of GX 304-1 . The in- 
creased solar activity during the observations led to the reduction 
in the exposure time and availability of the instruments as can be 
seen from the exposure columns of Table 1 . 

The INTEGRAL observations are indicated in Fig 1 , 
which shows the Swift/BAT light curve of GX 304—1 '. The 
INTEGRAL monitoring has an excellent coverage of the entire 
outburst, providing a rare opportunity to follow the evolution of 
the outburst from the early rising phase to the late decay phase. 

For our analysis, we used data from the ISGRI detector of 
IBIS, which is sensitive in the 20-300 keV energy range, JEM-X, 
and SPI. Standard data processing was performed with version 9 
of the Offline Science Analysis (OSA) software. We performed 
an additional gain correction of the ISGRI energy scale based 

1 Wc used the Swift/BAT transient monitor results provided by the 

Swift/BAT team 
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Fig. 1 . INTEGRAL observations of GX 304—1 (horizontal bars) 
superposed on the Swift/BAT light curve of the source during its 
outburst in January-February 2012 (-022 units of the vertical 
axis corresponds to one Crab). 

on the background spectral lines of Tungsten (the modified re- 
sponse files based on the nearest Crab observations were used). 

3. Spectral analysis 

Figure 1 shows that the X-ray flux of GX 304-1 during the ob- 
served part of the outburst changed by an order of magnitude. 
This allowed a detailed study of the luminosity dependence of 
the source’s broad-band X-ray spectrum. X-ray pulsations with 
a period of -274.9 s were detected in all INTEGRAL observa- 
tions. This value is roughly consistent with the known pulse pe- 
riod. A detailed timing analysis is not part of the present work 
and will be presented elsewhere. 

In all observations, the X-ray continuum could be closely 
modeled by a standard power-law/cutoff function (flux oc 
£T r exp [£/£fo!dL where E is the photon energy, T is the pho- 
ton index, and E(m is the exponential roll-off parameter) modi- 
fied by photo-electric absorption at low energies. In addition, the 
spectra showed a cyclotron resonant scattering feature around 
~50keV in absorption, which was modeled with a multiplica- 
tive absorption line with a Gaussian optical depth profile 1(E) = 
fcont (E) ■ e -0 ® , where G(E) = -r cyc /( V2rr<r cyc ) - exp[-0.5(E - 
Ecyc) 1 2 /o^y C ] , 4 om(E) is the continuum function, E^, cr eye, and 
Tcjc are the centroid energy, width, and optical depth of the line, 
respectively. The line is clearly detected in ISGRI and SPI data 
separately, as shown by the residual plots in Fig. 2. The inclusion 
of the absorption line in the model leads to an improvement in 
the reduced* 2 from 3.10 for 145 d.ojf. (with the corresponding 
null-hypothesis probability of only -10' 32 ) to 0.75 for 142 d.of. 
(null-hypothesis probability >0.9). The energy of the line is con- 
sistent with that reported by Yamamoto et al. (201 1) based on 
the RXTE and Suzaku observations. Wc also included an addi- 
tive Gaussian component to model the Fe K a fluorescence emis- 
sion line around 6.4 keV. The inclusion of the line reduces the 
residuals around 6 keV. The corresponding improvement of the 
reduced* 2 is, however, marginal: from, 0.80 (143 d.oi.) to 0.75 
(142 d.of.). The presence of the Fe line is, therefore, question- 
able. 

Both the continuum and cyclotron line parameters vary sys- 
tematically during the outburst. Here, we focus on the evolution 
of the cyclotron line energy E c yc, to establish a possible correla- 
tion of Eeye with flux, similar to that found by Yamamoto et al . 
(2011). To characterize this variability, we used only observa- 
tions where data from all three INTEGRAL X-ray instruments 
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Fig. 2. Example of the INTEGRAL spectrum of GX 304-1 (rev- 
olution 1132) with a simultaneous fit of data from JEM-X, 
ISGRI, and SPI with a power-law/cutoff model including a cy- 
clotron absorption line (see text) (a), the residuals for a fit with 
the model without the cyclotron line for JEM-X + ISGRI (b), 
ar.d JEM-X + SPI (c), and with the model where the line is in- 
cluded for data from all instruments (d). 



Fig. 3. Evolution of the cyclotron line centroid energy E cyc 
throughout the outburst of GX304— 1 as measured with 
INTEGRAL. The vertical error bars indicate lfr-uncertainties. 
The horizontal error bars indicate the time intervals of the obser- 
vations. The dotted line represents the re-scaled Swifl/BAT light 
curve. 


were available. This filtering allowed us to obtain a homoge- 
neous set of data and, hence minimize possible systematic ef- 
fects. Following this selection, data from the orbits 1133 and 
1135 were excluded. The JEM-X data collected during the orbit 
1134 (peak of the outburst) have been flagged as “bad” by the 
instrument team owing to the impossibility to provide a precise 
energy calibration. In our analysis, however, we use the JEM- 
X data to restrict the low-energy part of the broad-band X-ray 
continuum of the source, for which a precise energy scale is not 
very critical . Therefore, after consulting the JEM-X team, we re- 
introduced the JEM-X data from the revolution 1 134 in our anal- 
ysis. We checked, however, that inclusion of the JEM-X data of 
rev. 1 1 34 to the corresponding spectral fit does not significantly 
influence the measured energy of the cyclotron line (the focus of 
this work), but only affects its uncertainty. Table 2 summarizes 
the best-fit spectral parameters achieved in each observations. 


Lx/10 37 erg s 1 



Fig. 4. Cyclotron line centroid energy £ cyc as a function of the 
logarithm of flux in the 4-80 keV range. The error bars indicate 
1 (^-uncertainties (the flux uncertainties are smaller than the sym- 
bol size). The top X-axis represents the corresponding isotropic 
source luminosity assuming a distance of 2.4 kpc. The dotted 
line shows a linear fit to the £cy C -Iogjo(Flux) dependence. 


A clear systematic variation in the line energy over the out- 
burst is evident in Fig. 3, which shows that Ecy C generally fol- 
lows the X-ray flux. To assess the interdependence of the two 
parameters, we plotted Ecyc as a function of the X-ray flux in the 
4-80 keV range measured with INTEGRAL in the respective ob- 
servations (Fig. 4) . The plot shows a positive correlation between 
the two values. A linear fit to the dependence of E cyc on the loga- 
rithm of flux (dotted line in the plot) reveals a slope of 4.97±1 .12 
kcV/logio(erg cm -2 s -1 ). The standard linear correlation analysis 
of the Ecyc-logio(Flux) dependence yields a Pearson’s correla- 
tion coefficient of 0.88 with a probability of obtaining the corre- 
lation by chance of ~0.01 (one-sided). 

We verified whether our spectral fits contain artificial (model- 
driven) dependences of the cyclotron line energy E^ and other 
model parameters usings 2 -contour plots. No significant model- 
driven dependences between Ecy C and any of the continuum pa- 
rameters were found. The line energy was, however, found to be 
somewhat coupled to the line width cr c yc and its central optical 
depth T cyc . Nevertheless, the^-minima and the confidence inter- 
vals could be clearly identified and are separated for the different 
observations as shown in Fig. 5. The plot shows the contours for 
the /-eye/ tr eye pair. The EcycAcyc contours look similar. 

To check whether the Ecyc/flux correlation is related to in- 
strumental effects, we performed spectral fits using only the 
JEM-X and SPI data (excluding ISGRI). To verify whether the 
correlation depends on the choice of the spectral model, we fit 
the data using a Lorentzian line profile instead of a Gaussian 
one. We also tried two alternative continuum functions: XSPEC 
powerlawxhighecut and compTT models. In the former model, 
the highecut component controls the exponential roll-off. In ad- 
dition to Efoid, this component has an additional parameter - the 
cutoff energy E cut0 ff, above which the spectrum is affected by the 
roll-off. In our fits, E cuto a stays between a few and ~10keV. In 
all cases, the positive Ecyc/flux correlation was reproduced. We 
conclude, therefore, that the reported correlation arises from the 
source’s behavior and reflects real physics. 

4. Discussion and conclusions 

The positive correlation between the cyclotron line centroid en- 
ergy and the flux found with INTEGRAL confirms the claim 
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Table 2. Best-fit spectral parameters with the corresponding lcr-uncertainties for the INTEGRAL observations used in this work. 
Tie luminosity is provided assuming a distance of 2.4kpc and an isotropic emission diagram. 


Revolution 

1131 

1132 

1134 

1136 

1137 

1138 

r 


0.83 

0-93!““ 

I-HdS 

1 Tt+0.06 
0.08 

1 56*° 13 

** -005 

Efow, keV 


15-5^ 

'4.5^ 

16.9;$ 

19-5!!i 

28.6!™ 

Ecyc, keV 

52.98-* 

54.39!™ 

55.7»^ 

51.44^ 

51.70^ 

48.37! 

<r cjc. keV 

5 4540.62 

0.58 

7 - % !§3 

8-45^” 

5.60!^ 

5 44+0.36 

4.80!™} 

Tcyc 

9.4 +i | 


13.3!^ 

8.8!™ 

8.9!|* 

5.7!!I 

Flux/10' 8 ergs' 1 cm' 2 

0-375-™ 

0011 now 

l-630:“;°li 

0.729!!™ 

0.340!“;™ 

0 151 +0<K)2 
U - 1>:,, -0002 

Lx/lO 37 erg s ' 1 

0.26 

0.56 

1.12 

030 

0.23 

0.10 

xl,Jd.oS. 

1.1/164 

0.7/164 

13/164 

1.1/164 

0.9/164 

1.0/164 



Fig. 5. ^-contour plots of the parameter pair E^Jcr^ for a 
few selected observations. The contours correspond to x^ in + 1 .0 
(the projections of this contour to the parameter axes conespond 
to the 68%-uncertainty for one parameter of interest), x^ ia +2.3 
(68%-uncertainty for two parameters of interest), and +4.61 
(90%-uncertainty for two parameters of interest). The respective 
orbit numbers are indicated. 

of this correlation by Yamamoto etal. (2011). The revealed 
dependence establishes GX 304-1 as the third member of a 
slowly emerging class of accreting pulsars showing a positive 
Fcyc/flux correlation with the other members being HerX-1 
(Staubertetal. 2007) and possibly A 0535+26 (Klochkov etal. 
2011), for which the positive correlation has so far been only 
established in the pulse-to-puise analysis. The opposite (neg- 
ative) correlation between the line energy and the flux was 
found in V0332+53 and 4U0115+63 (e.g., Tsygankovctal. 
2007, 2010). According to discussions in Staubert et al. (2007), 
Klochkov et al. (201 1), and Becker et al. (2012), these two types 
of dependences reflect two different regimes of accretion. A 
particular regime is realized in a source depending on whether 
its X-ray luminosity L is above or below a critical luminos- 
ity It, which corresponds to the local Eddington luminosity at 
the X-ray emitting structure(s) on/above the neutron star sur- 
face. In accreting pulsars radiating above L? (“super-critical” 
sources), infalling matter is decelerated in a radiative shock, 
whose height is believed to increase with L, i.e., drift towards 
an area with a lower B-field strength. The opposite behav- 
ior probably occurs in sources radiating below or close to 
(‘'sub-critical” sources), where infalling matter is stopped by 
the Coulomb drag and collective plasma effects rather than in 


a radiative shock. As discussed in Staubertetal. (2007) and 
Becker etal. (2012), the height of the emitting region decreases 
with increasing luminosity owing to a corresponding increase 
in ram pressure of the infalling material, leading to a posi- 
tive E c yc/flux correlation. The critical luminosity Lc depends 
on the parameters of the accreting neutron star, but is gener- 
ally around a few times ~10 37 ergs" 1 (Basko & Sunyaev 1976; 
Staubert et al. 2007; Becker et al. 2012). Assuming a distance of 
2.4kpc (Parkes et al. 1980), the X-ray luminosity of GX 304-1 
in the 4-80 keV range during the reported INTEGRAL observa- 
tions varies between ~l.lxl0 36 ergs -1 and ~1 .13xlO ,7 ergs 1 . 
Thus, according to the described picture, the source should be- 
long to the class of “sub-critical” sources, for which a positive 
Ecyc/flux correlation is expected. The reported observations are, 
therefore, in agreement with the idea of two accretion regimes 
and increases the yet very small sample of accreting pulsars with 
established Ecyc/flux correlations. 
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